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基于 二 进 制 烟花 优化 算法 的 认 知 无 线 网 络 频谱 分 配 
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(贵州 大 学 大 数据 与 信息 工程 学 院 , 贵阳 550025) 


摘 要 : 针对 当前 无 线 频谱 资源 稀缺 和 利用 率 低 的 问题 ， 提 出 一 种 基于 二 进 制 烟花 优化 算法 的 频谱 分 配方 法 。 每 个 烟 
花 个 体 进行 分 布 式 爆炸 搜索 ， 并 对 最 优 烟花 的 爆炸 半径 采用 改进 公式 动态 更 新 ; 在 变异 环节 中 ， 针 对 粒子 间 信 息 交 流 
不 足 的 缺点 ， 引 入 遗传 算法 的 交叉 变异 算 子 ， 进 一 步 增强 种 群 多 样 性 ; 对 选 出 的 最 优 个 体 使 用 Metropolis 准则 进行 模 
拟 退 火 扰动 ,避免 陷入 局 部 最 优 。 仿 真实 验 表 明 ， 二 进 制 烟 花 优 化 算法 在 认 知 无 线 网 络 的 频谱 分 配 中 具有 寻 优 精度 高 ， 
收敛 速度 快 的 特点 ， 较 好 地 实现 网 络 效益 和 用 户 比例 公平 性 的 最 大 化 。 
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Spectrum allocation of cognitive radio network based on 
binary fireworks optimization algorithm 


Yang Jun, Zhang Damin', He Qing, Pan Zhiyuan 
(College of Big Data & Information Engineering Guizhou University, Guiyang 550025, China) 


Abstract: According to the scarcity of wireless spectrum resources and low utilization rate, this paper presented a spectrum 


allocation method based on binary fireworks optimization algorithm. Each firework individual carries out distributed explosion 


search, and the method proposed an improved formula to update the explosion radius of the optimal fireworks dynamically. 


For the lack of information exchange between particles in the variation process, the algorithm introduced the crossover 


| operator and mutation operator of the genetic algorithm to further enhance the population diversity. In order to avoid falling 
into a local optimum, the algorithm used simulated annealing with Metropolis criterion to perturb the selected optimal 
individuals. The simulation experiment shows that the method has the characteristics of high optimization accuracy and high 
convergence speed. It can achieve the maximization of the network efficiency and the proportional fairness among different 
users. 
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引用 源 。。 
前 ， 国 内 外 的 学 者 们 对 认 知 无 线 网 络 的 频谱 分 配 进 行 了 
于 无 线 通 信 技 术 的 迅速 发 展 ， 无 线 频谱 资源 稀缺 问题 日 广泛 的 研究 ,针对 认 知 无 线 网 络 的 体系 架构 和 频谱 接 入 等 技术 ， 
突出 ， 并 且 频 谱 利 用 率 低下 ， 存 在 较 多 的 频谱 空洞 * As 现 有 的 频谱 分 配方 法 主要 有 博弈 论 错误 ! 未 找到 引用 源 。， 竞价 
" 。 认 知 无 线 网 络 "" (cognitive radio network,，CRN) 是 一 种 具 拍卖 理论 第 :被 到 用 源 ' 和 图 论 着 色 理 论 种 洪 ! 和 引用 源 ' 等。 考虑 到 CRN 
有 认 知 过 程 的 网 络 ， 它 能 分 辨 当前 的 网 络 状态 ， 然 后 根据 这 些 节点 之 间 的 干扰 关系 ， 图 论 着 色 理 论 可 以 形象 地 表示 出 用 户 之 
状态 进行 规划 、 决 策 和 响应 。 针 对 认 知 无 线 网 络 提出 的 动态 频 间 的 地 理 位 置 ， 从 而 实现 频谱 资源 的 规划 与 管理 ， 是 一 种 方便 
谱 分 配方 法 可 以 很 好 的 解决 以 上 问题 ， 其 核心 思想 是 在 不 影响 有 效 的 方法 。 由 于 认 知 无 线 网 络 的 频谱 分 配 是 在 考虑 用 户 公平 
主 用 户 《 授 权 用户 ) 正常 通信 的 情况 下 ， 次 用 户 《 认 知 用 生 的 约束 条 件 下 进行 ， 是 一 种 典型 的 NP 难 问题 ， 所 以 也 可 使 
可 以 通过 实时 感知 技术 机 会 使 用 主 用 户 的 频谱 ， 从 而 有 效 的 分 和 群 智能 算法 用 于 解决 。 近 年 来 ， 学 者 们 将 图 论 着 色 理 论 与 智 
配 可 用 的 频谱 资源 , 使 得 网 络 整体 性 能 达到 最 优 错误 ! 未 找到 能 优化 算法 相 结 合 ， 并 应 用 到 频谱 分 配 当 中 ， 取 得 了 很 多 的 成 
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果 。 文 献 [6，7] 提 出 基于 遗传 (genetic algorithm，GA) 和 量子 
遗传 算法 的 频谱 分 配 解决 方案 ， 但 是 遗传 和 量子 遗传 算法 都 容 
易 早 熟 收 化 ， 影 响 系 统 的 最 终 寻 优 结果 。 文 献 [8，9] 在 认 知 无 
线 网 络 的 频谱 分 配 中 引入 了 免疫 克隆 选择 算法 (immune clone 
selection, ICS), 很 好 地 解决 了 种 群 进化 后 期 缺乏 多 样 性 的 问题 ， 
但 是 收敛 速度 不 够 快 。 文 献 [10] 采 用 粒子 群 算法 (particle swarm 
optimization，PSO) 进 行 频谱 分 配 ， 虽 然 有 较 快 的 收敛 速度 ， 但 


圆 形 区 域 .假定 用 户 间 的 干扰 仅 由 地 理 位 置 上 的 相互 距离 决定 。 
如 果 主 用 户 i 和 认 知 用 户 4 或 者 认 知 用 户 n 和 k(k=1,2,.…N) 在 
同一 频谱 m 上 的 覆盖 范围 存在 重 闭 则 表明 将 产生 干扰 。 
1.2 数学 建 模 
由 于 频谱 分 配 的 周期 相对 于 网 络 环境 的 变化 时 间 很 短 ， 所 
以 假定 在 进行 一 次 频谱 分 配 的 过 程 中 ， 网 络 拓扑 结构 不 发 生 改 
变 。 认 知 无 线 网 络 的 频谱 分 配 可 以 转化 为 数学 模型 进行 表达 。 
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是 容易 陷入 局 部 最 优 。 除 此 之 外 ， 还 有 很 多 其 他 的 经 典 智 能 算 
法 在 CRN 频谱 分 配 问 题 中 同样 也 得 到 了 广泛 的 应 用 ， 如 蚁 群 
算法 错误 ! 未 找到 引用 源 。、 蜂 群 算法 错误 ! 未 找到 引用 
源 。、 布 谷 岛 算法 错误 ! 未 找到 引用 源 。 等 。 但 是 这 些 方法 
都 很 难 在 收敛 速度 和 寻 优 能 力 之 间 取 得 平衡 ， 进 而 会 影响 最 终 
的 频谱 分 配 结果 以 及 网 络 性 能 。 
本 文 以 图 论 着 色 频 谱 分 配 模型 为 基础 ， 以 网 络 效益 以 及 用 
户 之 间 的 比例 公平 性 为 优化 目标 ， 对 传统 二 进 制 烟花 算法 和 * 虹 
ae 的 爆炸 半径 和 变异 算 子 进行 改进 ， 并 加 入 
Metropolis 准则 进行 择优 ， 提 出 了 一 种 基于 二 进 制 烟花 优化 算 
法 (binary fireworks optimization algorithm, BFOA) 的 频谱 分 配 
解决 方案 。 通 过 实验 仿真 比较 分 析 ， 证 明了 BFOA 算法 解决 该 
问题 的 有 效 性 和 优越 性 。 


1 ” 认 知 无 线 网 络 频谱 分 配 模 型 


1.1 频谱 分 配 图 论 模 型 

如 图 1 所 示 ， 图 论 着 色 模 型 以 无 向 图 的 形式 诠释 认 知 无 线 
网 络 的 拓扑 结构 ， 展 现 用 户 之 间 的 地 理 位 置 关系 以 及 互相 存在 
的 干扰 。 


图 1 认 知 无 线 网 络 拓扑 结构 
图 中 ， 主 用 户 括号 内 为 其 授权 频谱 ， 认 知 用 户 括号 内 为 其 
感知 到 的 可 用 频谱 ， ri 和 用 分别 为 各 自 的 信号 覆盖 半径 。 

段 设 在 一 片 X*Y 的 区 域 中 ， 存 在 有 7 个 主 用 户 和 N 个 认 知 用 
户 , 且 其 地 理 位 置 是 随机 分 布 的 , 区 域内 可 用 频谱 的 个 数 为 W 。 
主 用 户 使 用 固定 的 授权 频谱 ， 认 知 用 户 通过 认 知 无 线 网 络 的 频 
谱 感 知 技术 感知 当前 主 用 户 的 可 用 频谱 , 并 机 会 接 入 进行 使 用 ， 
在 满足 频谱 分 配 规则 的 前 提 下 ， 认 知 用 户 允 许 同时 使 用 多 个 频 
谱 。 当 主 用 户 i(i=42.…7) 分配 的 授权 频谱 为 m(m=1,2.…M) 时 ， 
则 主 用 户 i 在 频谱 m 上 的 覆盖 范围 是 以 i 为 中 心 ， 序 , 为 半径 的 
圆 形 区 域 。 若 认 知 用 户 n(n=4%…N) 机 会 接 入 了 频谱 mw ， 则 认 
知 用 户 在 频谱 sw 上 的 履 盖 范围 是 以 7 为 中 心 ，mm, 为 半径 的 


可 由 以 下 和 矩阵 表示 : 
a) 可 用 频谱 矩阵 ={4e{0,j,, ， 表 示 认 知 用 户 7 是否 
可 以 使 用 4,=1 则 表示 可 以 ， 相 反 ，4=0 则 表示 不 
可 以 。 

日 效益 矩阵 B= 他 |。 =(ms) | ,表示 认 知 用 户 使 用 频 
谱 m (14%=1) 时 所 得 的 效益 ， 且 效益 为 正 。 一 般 将 网 络 带宽 
作为 频谱 效益 。 

c) 干 扰 矩 阵 C={cwuwjcuw sfo1} ， ，， 表 示 不 同 的 认 知 用 户 
nn 和 在 使 用 同一 频谱 m 时 是 否 会 存在 干扰 。 c=1 则 表示 存 
在 干扰 ，c4 =0 则 表示 不 存在 。 
d) 无 干扰 分 配 和 矩阵 4={o osfoI}， ， ， 表 示 在 可 用 频谱 
和 矩阵 民 和 干扰 矩阵 c 都 满足 的 条 件 下 ,频谱 m 是 否 可 以 分 配给 
认 知 用 户 n。a=1 表 示 可 以 ，%=0 则 表示 不 可 以 。 无 干扰 
分 配 和 矩阵 须 满足 如 下 约束 条 件 : 


dim * pm = 0 crm =1 


频谱 。 


VISnkSN ,1<m<M (1) 


最 终 的 无 干扰 分 配 和 矩阵 和 效益 矩阵 ， 可 以 得 出 频谱 分 配 


完成 后 网 络 中 每 个 认 知 用 户 所 得 的 效益 。 即 
R, = 2 *b,, (2) 
网 络 获得 的 总 效益 即 为 所 有 认 知 用 户 获 得 的 效益 总 和 : 
F ,= NR, 
区 (3) 
=2 20 *b 
由 网 络 获得 的 总 效益 不 难得 出 其 平均 效益 : 
1 N 
Fean = 
1 N M (4) 
-3 
为 了 衡量 每 次 频谱 分 配 结果 的 公平 性 ， 即 网 络 中 各 认 知 用 
户 获 得 的 效益 是 否 均衡 ,可 以 用 网 络 的 比例 公平 性 函数 表示 千 & 


未 找到 引用 源 。 。 
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(5) 
-人 so 
通过 以 上 的 描述 和 分 析 可 以 看 出 ， 网 络 获得 的 效益 和 用 户 


间 的 比例 公平 性 可 以 用 来 反映 和 评价 认 知 无 线 网 络 频 谱 分 配 的 
结果 ， 即 可 以 直接 作为 算法 执行 过 程 中 的 目标 函数 。 
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2 BFOA 算法 设计 


传统 烟花 算法 由 爆炸 算 子 、 变 异 算 子 和 选择 策略 组 成 ， 且 
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Chin 
: 基于 二 进 制 烟花 优化 算法 的 认 知 是 瑟 网 线 闫 说/ 配 


杨 军 ， 
的 爆炸 半径 进行 改进 。 引 入 自 适 应 搜索 策略 ， 使 其 前 期 侧重 于 


常用 于 连续 变量 的 问题 中 。 为 适应 离散 化 的 频谱 分 配 应 用 ， 提 
出 一 种 全 新 的 二 进 制 烟花 优化 算法 。 在 变异 算 子 中 针对 粒子 信 
息 交 流 不 足 引入 了 经 典 遗 传 算法 的 交叉 变异 操作 ， 进 一 步 增强 
了 种 群 的 多 样 性 。 同 时 对 选 出 的 最 优 解 以 及 将 要 进行 下 一 轮 爆 
炸 的 最 优 烟 花 运用 模拟 退火 中 的 Metropolis 准则 进行 随机 扰 
动 ， 起 到 避免 算法 陷入 局 部 最 优 的 作用 ， 并 提高 了 收敛 速度 。 
此 外 ， 将 每 轮 从 代 的 最 优 烟 花 固定 爆炸 半径 改 为 动态 递 碱 ， 利 
于 自 适应 搜索 到 全 局 最 优 解 。 
2.1 爆炸 算 子 

种 群 中 的 每 个 个 体 属于 一 个 烟花 ， 烟 花 在 发 生 爆 炸 时 会 产 
生 大 量 的 火花 。 在 频谱 分 配 问 题 中 ， 适 应 度 值 较 大 (高 质量 ) 
的 烟花 爆炸 范围 较 小 ， 产 生 的 火花 较 多 。 相 反 ， 适 应 度 值 较 小 
〈 低 质量 ) 的 烟花 爆炸 范围 就 较 大 ， 产 生 的 火花 也 较 少 。 前 者 
具有 很 好 的 小 范围 “挖掘 能 力 ”, 后 者 则 具有 强大 的 大 范围 “ 探 
索 能 力 ”。 因 此 ， 种 群 中 的 每 个 烟花 根据 各 自 匹配 的 适应 度 值 
进行 爆炸 ， 最 终 可 以 使 得 种 群 在 全 局 搜索 能 力 和 局 部 搜索 能 力 
之 间 达 到 平衡 。 

每 个 烟花 的 爆炸 半径 
和 种 群 中 其 他 烟花 的 适应 度 值 共同 决定 ， 
应 的 爆炸 半径 和 火花 数量 分 别 为 


™t 


t 


和 产生 的 火花 数量 由 自身 的 适应 度 值 
烟花 i 发 生 爆 炸 时 对 


A =ceil| hv sh nf (NH)+e ee 
Pfr fe)+e 

S =ceil| S* ff te 人 
2 (f(%)- fn)+e 


其 中 : ceil 表示 向 上 取 整 ，A4 为 二 进 制 编码 维 数 ，$ 为 常数 ,用 
于 控制 产生 的 火花 总 数 。z 为 种 群 规模 ，/。 和 分 别 为 种 群 
中 最 大 和 最 小 的 适应 度 值 。。 为 机 器 精度 ， 取 ==10* 。 
当 第 i 个 烟花 发 生 爆 炸 时 ， 产 生 的 第 7 个 火花 表示 如 下 : 
x =O(x,P,r) (8) 
@(%,P,r) 为 二 进 制 转换 算 子 ， 表 示 从 角 标 集合 P={1,2.… 有 外 中 
随机 选 出 才 个 元 素 , 并 在 二 进 制 编码 的 字符 串 中 对 选 出 的 角 
标 对 应 的 字符 进行 0-1 替换 。 其 中 好 为 爆炸 算 子 步 长 ; 
m=ceil(A*rand+e) (9) 
表示 第 i 个 烟花 的 第 ;个 火花 半径 , 满足 0<m <4。 
上 服从 均匀 分 布 的 随机 数 。 
2.2 ”爆炸 半径 的 改进 
每 个 烟花 的 爆炸 半径 大 小 都 与 自身 的 质量 密切 相关 ， 质 量 
最 优 的 烟花 半径 最 小 ， 且 其 周围 的 粒子 相 比 其 他 烟花 也 最 接近 
全 局 最 优 。 所 以 需要 对 其 加 强 深度 的 “挖掘 ”， 才 能 提高 目标 
优化 的 精度 ， 因 此 质量 最 优 的 烟花 爆炸 产生 的 火花 数量 也 是 最 
多 的 。 但 是 通过 式 6) 计算 烟花 的 爆炸 半径 时 ， 会 容易 发 现 ， 
最 优 烟花 常常 会 由 于 半径 值 过 于 偏 小 而 难以 产生 出 最 多 的 火花 


错误 ! 未 找到 引用 源 。 


针对 这 一 问题 ， 本 文 在 Logistic 函数 的 基础 上 对 最 优 烟 花 


rand 是 [0,1] 


全 局 搜索 ， 而 后 期 侧重 于 局 部 搜索 。 
的 S 型 曲线 ， 能 够 较 好 地 反映 事物 
积分 简化 可 得 到 : 


Logistic 函数 是 一 种 常见 
了 界 增长 的 规律 ， 通 过 对 其 


i (10) 

变形 可 得 最 优 烟 花 半 径 公 式 : 

4(8)= 二 Er (11) 
人 
为 进化 开始 和 结束 时 设 定 的 最 优 烟花 半径 。 G 为 进化 总 代数 ， 


0 为 控制 参数 ，0<a<1。 

式 (11) 可 以 看 出 ， 随 着 进化 代数 的 不 断 增加 ， 最 优 烟 
花 的 半径 呈 非 线性 递减 ， 前 半 周 期 半径 较 大 ， 下 降 较 快 ， 后 半 
周期 递减 速率 越 来 越 缓慢 ， 且 逐渐 趋向 于 最 小 半径 值 。 这 样 有 
利于 算法 在 迭代 初期 进行 空间 内 的 广泛 搜索 ， 而 末期 能 够 在 全 
步 精确 的 局 部 搜索 。 


局 最 优 附 近 进行 进 
2.3 改进 变异 算 子 

传统 烟花 算法 的 变异 算 子 通常 采用 高 斯 变异 ， 但 高 斯 变异 
引导 个 体 跳出 局 部 最 优 解 的 能 力 较 弱 ， 且 一 般 应 用 于 实数 编码 
问题 。 此 外 ， 烟 花 算 法 属于 一 种 并 行 弥漫 式 爆炸 搜索 的 算法 ， 
在 寻 优 过 程 中 粒子 间 缺 乏 信息 交流 ， 多 样 性 有 待 进 一 步 提高 。 
本 文 引 入 遗传 算法 中 的 二 进 制 交叉 变异 算 子 对 其 进行 改进 ， 加 
强 了 粒子 间 的 信息 传递 和 资源 交互 ， 更 好 地 增强 了 种 群 的 多 样 


当 烟 花 完成 爆炸 后 ， 搜 索 空间 中 存在 烟花 和 爆炸 火花 两 种 


元 素 。 具 体 实现 过 程 如 下 : 
a) 用 目标 函数 评价 所 有 烟花 和 爆炸 火花 的 适应 度 值 并 排 
序 。 
b) 从 中 选择 合适 的 个 体 以 概率 p. 和 p, 执行 交叉 变异 操作 : 
dim = ceil(A* rand) if rand < p. (12) 
Xia if de randperm(A, dim) 
Cs b 而 otherwise (13) 
= | 一 % if rand 和 py 014) 
Wy otherwise 


其 中 : 和 芒 分 别 表示 交叉 和 变异 后 第 i 个 个 体 第 a 维 的 值 。 
randperm(A,dim) 表示 从 维 数 4 中 随机 选择 dim 个 不 重复 的 维 参 
与 交叉 。 

参与 交叉 变异 的 元 素 由 两 部 分 组 成 ， 一 部 分 是 适应 度 值 排 
名 靠 前 的 个 体 ， 另 一 部 分 则 在 剩 下 的 群体 中 随机 选 出 相同 数量 
的 个 体 。 这 样 更 有 利于 产生 优质 个 体 ， 同 时 也 没有 放弃 适应 度 
差 的 个 体 身 上 的 优秀 基因 ， 增 强 了 个 体 之 间 的 信息 共享 。 
2.4 选择 策略 
以 上 过 程 结 束 后 ， 需 要 对 当前 候选 者 集合 进行 一 定 规则 的 
筛选 便于 将 种 群 中 的 优秀 信息 传递 给 下 一 代 ， 继 续 作 为 烟花 进 
行 爆炸 。 其 中 候选 者 集合 包括 烟花 ， 爆 炸 火 花 和 变异 火花 。 首 
先 ， 采 用 精英 保留 机 制 ， 将 适应 度 最 好 的 烟花 (火花 ) 确定 性 
地 保留 下 来 ,剩余 个 体 将 使 用 轮 盘 赌 的 方式 从 候选 集合 中 选取 ， 
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个 体 被 选中 的 概率 为 


D(%) 
P(%)= 工 
六 2(5) (15) 
Dl)= Day) (16) 


其 中 : D(x) 为 个 体 * 与 其 他 候选 者 之 间 的 海 明 距 离 之 和 ，K 为 
参与 轮 盘 赌 选择 的 候选 者 数量 。 
2.5 ”Metropolis 接受 准则 

为 了 使 算法 不 易 陷入 局 部 最 优 ， 本 文 在 烟花 完成 爆炸 变异 
后 引入 了 模拟 退火 算法 的 Metropolis 准则 ， 对 选 出 的 精英 个 体 
进行 进一步 优化 。 模 拟 退 火 算 法 通过 模拟 热力 学 中 经 典 粒子 系 
统 的 降温 过 程 来 求解 组 合 优化 问题 中 的 极 值 问题 。 其 核心 
是 对 当前 的 状态 进行 随机 扰动 产生 一 个 新 状态 ， 并 通 


有 记 


注 漳 


Metropolis 准则 判断 是 否 接 受 该 新 状态 。 由 于 频谱 分 配 问题 属 
于 求解 最 大 值 问 题 ， 故 Metropolis 准则 如 下 : 
1 f(x)<f(%) 


p= -| ee (17) 
其 中 : 为 当前 解 ，% 为 新 产生 的 解 ，f (x) 为 目标 函数 ，7 为 
初始 温度 。Metropolis 准则 不 但 能 够 接受 优化 解 ， 而 且 还 能 够 
以 一 定 的 概率 接受 恶化 解 。 当 温度 较 高 时 ， 算 法 接受 恶化 解 的 
概率 较 大 ; 随 着 温度 的 递减 ， 接 受 恶化 解 的 概率 逐渐 减 小 ， 直 
至 最 后 趋 于 0。 从 而 更 有 利于 算法 跳出 局 部 极 值 以 趋向 于 全 局 
最 优 。 本 文采 用 文献 错误 ! 未 找到 引用 源 。 中 的 降温 函数 来 
制 温度 的 下 降幅 度 : 

TH+D=T(D*B (18) 
其 中 : 0<p6<1 为 降温 系数 ，4 为 当前 迭代 次 数 。 通 过 此 函数 进 
于 降温, 可 以 达到 退火 初期 降温 快 , 而 在 后 期 下 降 缓慢 的 效果 。 
且 ， 为 避免 在 退火 过 程 中 因 概 率 接受 恶化 解 而 错失 当前 为 止 
的 最 优 解 ， 在 迭代 过 程 中 添加 了 “best so far” 记 忆 单 


3 ”基于 BFOA 的 频谱 分 配方 案 


在 本 文 提出 的 基于 二 进 制 烟花 优化 算法 的 频谱 分 配方 法 
中 ， 每 个 烟花 和 产生 的 火花 (爆炸 火花 、 变 异 火花 ) 代表 一 种 
可 能 的 频谱 分 配 策 略 。 通 过 对 烟花 和 火花 的 干扰 约束 处 理 后 ， 


HY 


—、 


Sl 


o 


Xe, 


用 相应 的 目标 函数 进行 评价 从 而 最 终 选 出 最 优 的 频谱 分 配方 
案 。 
3.1 种 群 编码 


认 知 无 线 网 络 的 频谱 分 配 问题 采用 二 进 制 编码 的 方式 。 当 
可 用 频谱 矩阵 工 中 的 元 素 关 =1 时 ， 经 过 约束 条 件 的 处 理 ， 对 
应 的 无 干扰 分 配 和 矩阵 4 中 元 素 @ 可 能 为 0 也 可 能 为 1。 但 
bm=0 时 ，@w 必 为 0， 表 示 频 谱 从 刚 开 始 就 不 可 用 。 因 此 ， 可 
以 只 对 工 和 矩阵 中 为 1 的 元 素 进行 编码 ， 有 助 于 提高 计算 速度 。 
每 个 参与 计算 的 烟花 和 火花 个 体 其 实 是 频谱 分 配 矩 阵 到 二 进 制 
字符 串 的 映射 ， 是 一 种 可 能 的 频谱 分 配方 案 。 字 符 串 长 度 即 为 


~ 


- 关 >4。， 编 码 顺序 遵循 先 " 递增 ， 后 mm 递增 的 规则 。 
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杨 军 ， 等 : 基于 二 进 制 烟花 优化 算法 的 认 知 无 线 网 络 频 谱 分 配 
3.2 干扰 约束 处 理 
在 频谱 分 配 过 程 中 ， 并 非 每 一 种 分 配方 案 都 是 可 行 的 。 因 
为 可 能 会 出 现 多 个 认 知 用 户 使 用 同一 频谱 时 产生 干扰 冲突 的 问 
题 ， 故 需要 进行 干扰 约束 处 理工 作 。 当 认 知 用 户 n 和 在 使 用 
频谱 m 时 存在 干扰 ， 即 干扰 矩阵 元 素 Gxm=1 ， 则 频谱 m 只 允许 
个 用 户 接 入 。 此 时 需 检查 分 配 和 矩阵 4 的 第 m 列 第 n 行 和 第 x 
行 元 素 ， 若 都 为 1 则 需 处 理 : 比较 两 行 的 元 素 之 和 sum(n) 和 
将 最 大 者 对 应 的 m 列 元素 置 换 为 0。 即 若 
sum(n)>sum(k) ， 则 将 第 4 行 m 列 改 为 0; 若 sun(n)=sum(k)， 
则 随机 选 一 改 为 0。 经 过 干扰 处 理 后 的 分 配 和 矩阵 4 就 是 满足 矩 
阵 L 和 Cc 约束 条 件 的 有 效 解 。 
3.3 ”算法 实现 流程 
基于 二 进 制 烟 花 优 化 算法 的 频谱 分 配 具体 实现 流程 如 下 : 

a) 生 成 矩阵 ,。 由 随机 生成 的 CRN 网 络 拓扑 结构 图 得 出 可 用 
频谱 矩阵 LK ， 效 益 矩 阵 B 以 及 干扰 矩阵 C 。 

b) 初 始 化 参数 。 输 入 种 群 规模 Z ， 最 大 友 代 次 数 G ， 交 又 
变异 概率 p. 和 p, ， 初 始 温 度 7 等 ， 并 随机 产生 初始 烟花 种 群 
小 =[ 芭 , 码 ,…, 蔚 ] ， 使 用 相应 的 目标 函数 计算 初始 种 群 的 适应 度 
值 。 


Spe 


sum(k) ， 


c) 终 止 条 件 判断 。 检 测 是 否 达到 最 大 友 代 次 数 G ， 若 达到 
则 退出 循环 ， 将 当前 最 优 烟 花 进行 分 配 矩 阵 映 射 ， 输 出 最 优 频 
谱 分 配方 案 ; 若 没有 则 继续 执行 d)。 

d) 爆 炸 操 作 。 根 据 计 算得 出 的 爆炸 半径 4 、4(8) 和 火花 
数量 5; 对 每 个 烟花 进行 爆炸 ， 生 成 对 应 的 火花 马 。 

e@) 变 异 操 作 。 从 烟花 和 产生 的 火花 集合 中 选 出 适应 度 值 排 
名 靠 前 的 50% Z 个 元 素 , 再 从 剩 下 的 个 体 中 随机 选 出 50% Z 个 
元 素 组 成 参与 变异 的 群体 , 以 交叉 概率 p. 和 变异 概率 p, 执行 交 
又 变异 操作 ， 生 成 变异 火花 。 
fj) 约束 处 理 与 择优 。 对 烟花 和 产生 的 爆炸 火花 、 变 异 火 花 
群体 进行 干扰 检查 和 约束 处 理 。 并 计算 适应 度 , 选 出 最 优 个 体 。 
g) 搜 索 全 局 最 优 。 利 用 Metropolis 接收 准则 对 选 出 的 最 优 
个 体 进行 扰动 更 新 并 约束 处 理 ， 进 一 步 搜索 将 进行 下 一 轮 爆炸 
的 全 局 最 优 烟 花 。 

h) 选 择 更 新 。 将 g) 中 选 出 的 最 优 烟花 对 了 f) 中 的 最 优 个 体 进 
行 蔡 换 ， 并 与 其 剩 下 个 体 组 成 集合 。 使 用 轮 盘 赌 的 方式 产生 其 
他 NW-1 个 进行 下 一 次 爆炸 的 烟花 种 群 。 

is=s+1， 转 至 c)。 


4 ”仿真 及 结果 分 析 


利用 MATLAB 建立 仿真 实验 环境 。 在 相同 的 网 络 拓扑 结 
构 基 础 上 ， 将 提出 的 二 进 制 烟 花 优化 算法 (BFOA) 与 传统 解 
决 该 问题 的 遗传 (GA)、 粒 子 群 (PSO) 和 免疫 克隆 选择 (ICS ) 
算法 进行 比较 ， 分 别 通过 系统 网 络 效益 、 用 户 的 比例 公平 性 来 
分 析 四 种 算法 的 性 能 各 下 并 到 引用 源 ，- 错 未 入 到 引用 汤 。。 以 及 当 网 络 环境 
中 关键 参数 发 生变 化 时 对 算法 性 能 产生 的 影响 。 验 证 BFOA 用 
于 解决 认 知 无 线 网 络 频谱 分 配 问题 的 可 行 性 与 优越 性 。 在 给 定 
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的 矩形 网 络 区 域 中 ， 


随机 分 布 / 


个 可 用 频谱 。 认 知 用 


Z =10， 烟 花 总 数 为 


100， ms =0.257， 


交叉 变异 概率 产 =0.6， p, =0.1， 
4.1 算法 优化 性 能 比较 


本 组 实验 将 比较 四 种 算法 的 寻 优 能 力 和 收敛 性 能 。 图 2 
户 、 认 知 用 户 ) 位 置 分 布 图 

进行 频谱 分 配 时 网 络 总 效益 的 寻 优 结果 。 主 用 户 1、 认 知 用 
量 均 设置 为 15。 表 1 对 40 次 的 寻 优 结果 分 
行 求 和 和 均值 比较 。 从 图 2 和 表 1 中 可 看 出 ， 本 文 提 出 的 
获得 的 网 络 总 效益 最 大 ， 且 与 其 他 算法 的 


各 算法 在 40 次 不 同 


N 、 可 用 频谱 W 数 


演 


二 进 制 烟 花 优化 算法 
效益 差 值 较 明 显 。 


户 通过 动态 接 入 与 3 
量 实验 得 出 如 下 的 参数 设置 : 进化 总 代数 G =200， 种 群 规模 
rn =0.01 1 ， & =0.05， 


的 用 户主 用 


个 主 用 户 ，N 个 认 知 用 户 和 4 


E 用 户 共 享 频谱 。 经 过 


降温 系数 8 =0.98 。 


3 了 部 


图 2 不 同 用 户 位 置 分 布 图 下 的 网 络 总 效益 
表 1 40 次 网 络 总 效益 比较 
算法 40 次 效益 总 和 40 次 效益 均值 
BFOA 15554.4 388.86 
PSO 14972.0 374.30 
ICS 14398.9 359.97 
GA 14138.6 353.47 


图 3 和 4 分 别 是 四 种 算法 
益 和 用 户 比 例 公 平 性 的 优化 过 程 。 主 


在 进行 


i 


频谱 M 数量 均 设 置 为 10。1 


图 可 知 ， 


次 频谱 分 配 时 网 络 总 效 


主 用 户 1、 认 知 用 户 N、 可 


BFOA 算法 在 优化 过 程 


精度 上 看 ，BFOA 算法 


中 整体 性 能 高 于 其 他 三 种 算法 。 


从 寻 优 


寻 得 的 最 优 值 最 大 ， 
最 终 的 网 络 总 效益 和 比例 公平 
分 别 高 出 8.3% 和 8.9%。 由 于 


法 ， 具 有 执行 效率 
经 高 出 其 他 算法 。 


其 次 分 别 


是 PSO、 
性 收敛 值 


高 的 特点 ， 


4 中 分 别 在 80 代 和 60 代 左 右 


度 ， 且 能 很 好 地 跳 


出 局 部 最 优 。 


烟花 算法 
所 以 在 进 


ICS、GA 算法 。BFOA 
比 效 果 表 现 最 差 的 GA 
属于 分 布 式 爆炸 搜索 算 
化 开始 时 的 寻 优 值 就 已 


从 收敛 性 能 上 来 看 ，BFOA 算法 在 图 3 和 图 


就 已 经 收敛 ， 具 有 最 快 的 收敛 速 


这 是 天 


选 出 的 最 优 烟花 进 
快 地 寻找 到 全 局 最 


行 了 一 定 程 
优 解 


度 的 模拟 退 


为 在 每 次 迭代 过 程 中 对 
火 扰动 ， 使 其 能 够 更 


杨 


双全 全 测 。 


Chi 
军 ， 等 : 基于 二 VE 


0 20 40 60 80 100 120 140 160 180 200 
迭代 次 数 g 


3 网 络 总 效益 优化 过 程 


4.2 ”可 用 频 训 
当 网 络 环境 中 的 可 用 


20 40 60 80 10 120 140 160 180 200 
选 代 次 数 g 


四 4 比例 公平 性 优化 过 各 
普 数量 对 算法 影响 的 比较 


频谱 数量 发 生变 化 时 ， 算 法 的 性 能 将 


会 受到 相应 的 影响 。 图 5 和 6 分 别 验证 了 在 不 同 频谱 数 时 网 络 
平均 效益 和 比例 公平 性 的 变化 情况 。 这 里 认 知 用 户 数 取 N =15， 


主 用 户 了 I 逢 


[可 用 频谱 W 数量 变化 区 间 为 5 到 30。 在 图 5、6 中 ， 


随 着 可 用 频谱 数量 的 增加 ， 系 统 网 络 平均 效益 和 比例 公平 性 都 


在 不 断 提升 。 这 是 由 于 当 频 谱 数 量 增加 时 ， 


法 得 出 的 | 


比例 公平 


有 效 地 证 明了 BFOA 算法 的 优越 性 。 


干扰 都 将 会 


明显 的 看 出 , 本 文 的 BFOA 算法 条 


改 益 差 距 不 是 很 大 ， 后 


生 远大 于 ICS 和 GA 算法 , 有 BFO 


r 


15 20 
可 用 频谱 数 M 


系统 的 频谱 冲突 和 


减轻 。 在 变化 的 过 程 中 ， 前 期 频谱 数 较 少 ， 四 种 算 
期 频谱 数 越 来 越 多 时 ， 可 以 较 
1 PSO 算法 获得 的 平均 效益 和 


A 的 值 比 PSO 更 大 ， 


5 可 用 频谱 数量 对 网 络 平均 效益 的 影响 
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4.3 


为 了 进 
上 ， 提 出 当 
会 产生 怎样 


比例 公平 性 


15 20 
可 用 频谱 数 M 


图 6 可 


频谱 数量 对 比例 公平 
认 知 用 户 数量 对 算法 影 
步 对 比 四 种 算法 的 性 能 ， 在 以 上 比较 分 析 的 基础 
区 域 中 的 认 知 用 户 数量 发 生 改 变 时 ， 对 相关 算法 将 


下 


性 


的 影响 


响 的 比较 


名 


的 影响 。 


到 30 


生 随 认 知 用 
户 了 和 可 用 频谱 数 M 均 取 15， 
网 络 环境 中 


户 数 变化 而 


妈 7 和 8 
变化 的 


展示 了 网 络 平均 效益 和 比例 公平 


认 知 


线 。 在 本 次 对 比 实验 中 ， 主 用 


j 户 数 N 的 变化 区 间 为 5 


Ne 


可 知 ， 当 


o 


随 着 认 知 用 
之 间 的 干扰 冲突 也 随 之 增加 
在 不 断 衰减 。 但 BFOA 算法 


户 数量 的 上 升 ， 


PF 的 主 


网 络 环境 中 的 负荷 加 重 ， 认 知 用 户 


户 和 可 用 频谱 数 固定 ， 


， 系 统 的 平均 效益 和 比例 公平 性 都 
改 获 的 效益 值 和 公平 性 仍然 是 四 种 


算法 中 最 高 的 ， 与 PSO、ICS、GA 相 比 更 能 实现 认 知 无 线 网 络 
0 


的 最 优 频 谱 分 配 ， 目 


5 


在 认 久 


用 户 数量 较 多 时 更 明显 。 


网 络 平均 效益 
二 全 


认 知 用 户 数 N 


户 数量 对 网 络 习 


F 均 效 


益 的 影响 


认 则 用 户 数 N 
图 8 认 知 用 户 数量 对 比例 公平 性 
结束 语 


的 影响 


针对 当今 无 线 通信 技术 中 频谱 资源 稀缺 以 及 利用 率 不 高 的 


问题 ， 


一 中 


本 文 将 二 进 


络 的 频谱 分 配 问题 相 结 合 ， 


后 的 


陷入 


得 
待 蝇 


判 烟花 算法 进行 优化 改进 ， 并 与 认 知 无 线 网 
了 较 好 的 频谱 分 配 结果 。 
进 制 烟花 算法 解决 了 传统 烟花 算法 中 收敛 速度 较 慢 ， 易 
局 部 最 优 解 的 问题 。 仿 真实 验 的 结果 分 析 表 明 ， 本 文 提 出 


优化 
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的 BFOA 算法 在 进行 频谱 分 配 时 , 与 其 他 三 种 算法 相 比 能 够 获 
得 更 好 的 网 络 效益 和 用 户 比例 公平 性 ， 从 而 提升 了 网 络 的 整体 
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